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Abstract 

We derive a Melvin Universe type solution describing a magnetic field permeating the whole 
Universe in gravity minimally coupled to any non-linear electromagnetic theory, including Born- 
Infcld Theory. For a large set of non-linear electrodynamics theories, our solution is complete 
and non-singular, as long as the magnetic field is sub-critical. We examine some properties of the 
solution; in particular there is a shift of the symmetry axis and a non-standard period along the 
orbits of the U(l) symmetry to avoid a conical singularity. We show these are consistent with 
the usual Dirac quantization condition for the magnetic flux. We find exact solutions describing 
propagation of waves in the 'generalized Melvin Universe' along the principal null directions of 
the electromagnetic field, where the Boillat and Einstein light-cone touch. By electric-magnetic 
duality we show that similar Melvin electric and dyonic universes can be obtained. 



1 Introduction 

The 'Melvin magnetic universe' is a simple solution to Einstein-Maxwell theory (or the bosonic sec- 
tor of N=2 Supergravity) describing a bundle of magnetic flux lines in gravitational-magnetostatic 
equilibrium (!]][]. Subsequent work has generalized the solution in different directions, including 
rotating and time dependent magnetic universes ||, black hole spacetimes which are asymptoti- 
cally Melvin 0], gravitational waves travelling in the magnetic universe || and a Melvin universe 
with a non-trivial axion and dilaton in heterotic string theory Q. Particularly interesting is the 

'Email address: gwgl@damtp.cam.ac.uk 
^Email address: car26@damtp.cam.ac.uk 

1 See also comments in [^| referring to work previous to Melvin's on the same solution. 
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Melvin Universe solution in Kaluza-Klein theory since it derives from a flat spacetime by perform- 
ing U(l) Kaluza-Klein reduction with a twist in the identifications pCj]. Recently [18| a version of 
the Kaluza-Klein Melvin Universe was derived from M-theory as a solution to either type IIA or 
type OA string theory. 

One of the most interesting applications of the Melvin Universe is the study of black hole pair 
creation in a magnetic field. Although this spacetime configuration is stable against arbitrarily 
large radial perturbations |?J , it is unstable against the non-perturbative process of black hole pair 
creation. This process was first discussed in based on the 'asymptotically' Melvin 'C-metric' 
type solution found by Ernst ||. The Euclideanized version of this metric is in fact the instanton 
describing the creation in the Melvin Universe of a pair of magnetically charged extreme Reissner- 
Nordstrom black holes with opposite charges p0[ . In the Kaluza-Klein cases [21|, a second decay 
mode is by bubble nucleation. Moreover, and quite strikingly, the Ernst instanton is nothing else 
than the Kerr solution in the higher dimension [21]. In M-theory the instability produces pairs of 
D6-branes and anti-branes in the type IIA case, whereas for the type OA the spacetime decays via 
bubble nucleation. 

In this note we generalise the Melvin spacetime in yet a different direction - as a solution to 
gravity minimally coupled to any theory of non-linear electrodynamics (NLE) in four spacetime 
dimensions. Our approach relies on a limiting procedure in Einstein-Maxwell theory which trans- 
forms the Reissner-Nordstrom solution into the Melvin solution. Since the Reissner-Nordstrdm 



solution is known for an arbitrary theory of non- linear electrodynamics (see [11] which includes the 
original references), it is straightforward to generalise the procedure to Einstein-NLE and hence 
find the generalized Melvin solution. 

Our motivation stems both from the fact that it might help us understand phenomena like black 
hole pair creation in gravity coupled to NLE, most notably Born-Infeld theory, and in connection to 
our recent work on causality in string theory and NLE [12|. Therein we showed that the open string 
metric that has recently arisen in the work of Seiberg and Witten [ 13 1 is just the Boillat metric [|l4|] 
describing the propagation of fluctuation of the electromagnetic field in a non-trivial background 
electromagnetic field in NLE. The Boillat metric is not, in general, conformal to the Einstein or 
closed string metric. Hence, we expect light and gravitons to travel at different speeds, except along 
the principal null directions of the background electromagnetic field [^] . The generalized Melvin 
solution we obtain in this paper will allow us to give examples of gravitational waves travelling in 
a curved spacetime along the direction in which the Boillat and Einstein light-cone coincide. 

This paper is organised as follows. In section 2 we describe the formalism of non- linear electro- 
dynamics coupled to gravity and the spherically symmetric charged solution which will be essential 
for our derivation. In section 3 we describe the limiting procedure and derive the solution. In 
section 4 we discuss some properties of the solution and compare it to the Maxwell case. In par- 
ticular, we discuss the existence of curvature singularities beyond the critical field in Born-Infeld 
theory. A general shift of the symmetry axis is found. In order to avoid conical singularities a 
non-standard period must be imposed along the orbits of the azimuthal Killing vector field. We 
also discuss quantization of magnetic flux, showing that the standard condition is still obtained. In 
section 5 we find exact solutions describing gravitational waves in the generalized Melvin Universe. 
In section 6 we use electro-magnetic duality to derive electric and dyonic Melvin Universes. We 
close with a discussion. 



2 



2 General Relativity coupled to non-linear electrodynamics 



2.1 General Framework 

If we consider gravity minimally coupled to non-linear electrodynamics, the action is 

S= — 1 f d A x^R+^- [ d 4 x^g~L(x,y). (2.1) 



16ttG n J v u 4vr . 

The quantities x = -^F^F^ and y = \F^F^ V are the two only independent relativistic invariants 
constructed from the Maxwell field in four spacetime dimensions, and L(x,y) an arbitrary function 
of those variables. We will take Newton's constant Gn = 1- 

The equations of motion derived from this action are more easily written in terms of the Legendre 
dual description of non-linear electrodynamics (see, for instance, J12||). This involves introducing 
the tensor P^ defined by 

dL = -Ip^dF^. (2.2) 

P^ u coincides with F^u for Maxwell's theory. In general it reads 

Pfj,v = —{LxF^v + Ly-k Fpv), (2.3) 

where subscripts on L denote differentiation. The components of P^ are just D and H, the electric 
induction and magnetic field respectively and hence ( |2.3j ) are faced as the constitutive relations for 
the non-linear electromagnetic theory. We then denote by s and t the two independent invariants in 
terms of the dual variables, i.e., s = —\P^i/P^ u and t = —\P^ V * P^ v ■ The covariant Hamiltonian 
L, 

L = -lp^F^-L, (2.4) 

is a function of (s,t). In an orthonormal frame, L coincides with the usual Hamiltonian for purely 
electric configurations and minus the usual Hamiltonian for purely magnetic ones. The equations 
of motion for the coupled system are fll7| 



D M P^ = 0, iV " \ R 9»v = 8ttT^, (2.5) 



with the energy momentum tensor^] 

T ^ = ^ (LsP^9 af3 Pup + g^u (2sL s + tL t -L^y (2.6) 

2.2 Spherically symmetric solution 

The spherically symmetric solution has been known since the 1930's. One finds that Birkhoff's 
theorem holds. For a general theory, the electrically charged solution may be cast in the form (see 
11] which includes the original references) 



Einstein = " ( 1 " ^P) *? + ( 1 - ^P) ^ dr* + r\dP + sin 2 



p Qdt A dr 



(2.7) 



We use a mostly plus signature. The Maxwell Lagrangian is L = —x = s. 



3 



The latter is the radial component of the electric induction D. The only non-trivial invariant is 
then 

a = Q 2 /(2r 4 ). (2.8) 
The metric function that defines the solution obeys the differential equation 

dm(r) 



dr 



r z L(s), (2.9) 



for s as in fl2.8| ). These equations permit one to read off the behaviour of the solutions given the 
function L(s). But it is actually easier to implement another program: choose m(r) such that 
the geometry is non-singular and find the 'theory' that gives rise to such non-singular solution. 
This program has been implemented in some recent papers [16]. Therein, possibilities for m(r) are 



considered so as to approach the usual Reissner-Nordstrom solution at large r. Furthermore, the 
geometry is chosen as to yield no curvature singularities; all scalar polynomials in the curvature are 
bounded for small r. One then finds the theory for which such geometries are solution via (|2.9| ). 
Although one may construct interesting solutions in this fashion, resembling Bardeen tunnels, they 
derive from obscure theories of non-linear electrodynamics. 



3 The Melvin Universe 

We now derive the generalized Melvin Universe solution to non-linear electrodynamics coupled to 
gravity. We start by describing a limiting procedure that generates the Melvin Universe of Einstein- 
Maxwell theory from the Reissner-Nordstrom solution. Then we apply the same procedure to the 
non-linear electrodynamics case, and verify the solution by solving the Einstein equations. The 
fields we obtain take the form: 



d$2 = p 2 + ^Hp) ^ + j?Mpf dp2 + A(p)2 { _ dT2 + dZ% 



p 2 + ^Kp) 

Hp 



HpY 



where the functions h(p) and A(p) are given by 



Hp) 



_i 

S 4 



L(s')s'—ids' 



• A dp, 



_H2 

2A(p)-J 



Hp) = i + ^V- 



(3.1) 



(3.2) 



For Maxwell's theory of electromagnetism, L = s and so one may easily check that h(p) = 0. The 
expression for h{p) might seem ill-defined, since it involves quartic roots that should be evaluated 
at a negative quantity. It should be understood, however, that first one should evaluate the integral 
and only then replace s by its value. In fact the expression is just a restatement of the infinite sum 
( [3.131) . The constant H in the solution defines the magnetic field strength at the origin, which is, 
in general, no longer at p = as shown below. 
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3.1 Melvin Universe in Einstein-Maxwell theory 

Let us describe the limiting procedure that generates the Melvin universe starting with the Reissner- 
Nordstrom (RN) solution. Start by writing the RN solution using stereographic coordinates for the 
S 2 sections: 



, 9 i , 4r 2 2M Q 2 

(l + cc) r r 

Q 

F = --\dr A dt, 



(3.3) 



Perform the diffeomorphism given by 

C = p r = rA, t=|, (3.4) 

to get the RN solution in (t, f, £, £) coordinates. Next we want to take the following limit: 

A -» 0, M, Q oo, MA 3 = /ixed = m, QX 2 = fixed = q. (3.5) 

That the resulting spacetime configuration remains a solution requires proof. We will come back to 
this point below. Notice that in the limit the S 2 is 'flattened' into R 2 . This changes the symmetry 
from spherical to axial. Denote £ = x + iz and perform the final coordinate transformations and 
relabellings: 

2 T Z B 2 2 

t==%< Pjj2> x = z "2' Z= 2' r = 1 + X P = 

(3.6) 

q = iB, m = -B 2 /2. 

Notice the change from timelike to spacelike character of t. This will change the electromagnetic 
field from electric to magnetic. The resulting configuration is the Melvin Universe ||]: 

ds 2 = P 2 k{ P )- 2 d4> 2 + A{p) 2 [-(IT 2 + dZ 2 + dp 2 ] , 

(3.7) 

F = BpA(p)- 2 dp A dep. 

The constant B defining the solution is the magnetic induction strength at the origin, since B = 
BA(p)- 2 d z . 

3.2 Generalized Melvin Universe 
3.2.1 The Limiting Procedure 

Now consider gravity coupled to general non-linear electrodynamics as described by ( |2.1| ). The 
generalized RN solution is (2.7). It takes the form fl3.3| ) with 



2M Itr' 1 ) , N 

V = l + 3.8 
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and F — ► P the dual Maxwell tensor. Equation (2.9) translates into 



,dl(r 



-1\ 



—L +r - l I(r- l )=2i*L( 

-1\ V } l 2r 4 



d{r~ x ) 



(3.9) 



for the function J(r _1 ). Next perform ( |3.4| ). In order to solve the differential equation we expand 
I(Af _1 ) in a power series (around spatial infinity) and L(s) in its Taylor series around zero (we are 
assuming L(s) to be analytic; we will be able to relax this assumption in the final result): 



KXr- 1 ) = Q 2 J2 /nCAf 1 )", Us) = J2 



n=0 



m=l 



ml 



d m L{s) 
ds m 



(3.10) 



s=0 



Requiring that L approach the value for Maxwell's theory in the weak field limit yields [dL/ds] s= o = 
1 and for the usual RN solution to hold one must have fo = 1. Substituting in fl3.9| ) and matching 
powers of Ar _1 we find 



oo 1 
n=0 v ' y 



Q 2 A 4 

(1 + n)!(l + 4n) V 2f= 4 
Taking the limit ( |3.5| ), which for the function / yields 



d n+1 L(s) 
ds n+1 



(3.11) 



lim A 4 /(Af-^ 

A->0,Q^oo,A 2 Q=g 



^ i 



n=0 



(1 + n)!(l + 4n) V2f 4 



,2 \ n 



9" +1 L(g) 



(3.12) 



s=0 



Finally perform (|3.6| ) with relabelling B ^ H to obtain the generalized Melvin magnetic universe 
(|3.1|). The function h(p) is 



-f I (- 

A,(l + n)!(l + 4n) V 



9 n+1 L(s) 



(3.13) 



s=0 



2A(p) 4 

It can also be expressed in the more useful integral form ( |3.2[ ). 
3.2.2 Solving the Einstein equations 

As we mentioned above, it is not clear that after performing the limit ( |3.5D one ends up with a 
solution of the equations of motion. The Maxwell equation in ( |2 .5[ ) is clearly obeyed by Q3.1| ) since it 
is completely identical to the usual Melvin solution and the determinant of the generalized Melvin 
metric is unchanged from the usual Melvin solution. So we just need to show the Einstein equations 
are solved. The Einstein tensor for the geometry ( p.l[ ) with A(p) given by ( |3.2j ) but unspecified 
h{p) is diagonal and so is the energy-momentum tensor. The pp and 4>4> equations of motion are 
equivalent, yielding the first order ordinary differential equation (O.D.E.) 



4s- 



, dh(s) 
ds 



s + sh(s) = —L(s). 



We are using the form of the only non- vanishing invariant 

H 2 



(3.14) 



(3.15) 
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to express h(p) as a function of s. It is easy to verify that h{p) given by ([^) (without specifying the 
value of s in terms of p) is the solution. Notice that this equation does not define the lower limit in 
the integral (^2); this is set by the requirement that one should recover the usual Melvin spacetime 
when L to correspond to Maxwell's theory. We can therefore drop the analyticity assumption for 
L(s). The tt and zz equations are also equivalent, yielding the second order O.D.E. 

- 8. 3 ^ - 14 S 2 ^ + ,(1 ~h) = 2sL s - L. (3.16) 
as* as 

This equation is not independent of the first order O.D.E. ( |3.14| ). Using the latter to find an expres- 
sion for L s and replacing this in the second order O.D.E. one recovers the first order equation. This 
is similar to what happens for solution (^]7]). Equation ( |2.9| ) comes from the tt or rr components 
of the Einstein equations, whereas the 69 and (fxj) components give a second order equation which 



is obtained by differentiating (2J}). Hence we have just shown that indeed the fields (|3.l| ) are a 
solution to the equations (|2.5[) . 



4 Properties of the solution and special cases 



We now discuss some properties of the solution ( |3.1| ). We will enphasize two particular types of 
theories of non-linear electrodynamics: 

i) The family described by the Lagrangian: 



L 



a n b 2 
-x + -x + -y , 



which has a dual Lagrangian (up to terms of higher order in the constants a, b), 



~ a 9 b 9 

L = s - -s 2 - -t 2 . 
2 2 



(4.1) 



(4.2) 



The Euler-Heisenberg theory (see for instance |L9|]) is the special case with a = 8a 2 /(45m 4 ) 
and b = 14a 2 /(45m 4 ), where m is the electron mass and a the fine structure constant. The 
function h(p) in (|3.2| ) becomes 

KP) = -^TTTI- (4-3) 



20A(p) 4 



ii) The Born-Infeld Lagrangian: 



L = -p(l- Vl + 2/3 2 x-/?V) , 



(4.4) 



which has an exact dual Lagrangian of the same form 

1 



L = p2 ( V 1 + 2/3 2 s - P 4 t 2 - 1 



(4.5) 

In string theory f3 = 2-Ka'. The integral ( |3.2| ) does not have a closed form and the best one 
can do in this case is to re-express h(p) in terms of the infinite sum ( |3.13| ): 

Hp) 



V i ( H 2 p 2 \ 1 

^(l + n)!(l + 4n) \2A( P r J 



H(2m-1) 



,m=l 



(4.6) 



Notice that both (|4.3|) for the Euler-Heisenberg case and (|4.6D are negative for all values of p. 
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4.1 Singularities and asymptotics 

The Einstein-Maxwell Melvin Universe ([T?]) represents a parallel bundle of magnetic flux held 
together by its own gravitational attraction. It is a totally non-singular configuration. It is easy to 
check the absence of singularities in scalar polynomials of the curvature: 

R = 0, R^IT = AB 4 A{p)- 8 , R^apR^ = 4B 4 A(p)~ 8 (3B 4: p 4: - 24B 2 p 2 + 80). (4.7) 

That R = follows from the conformal invariance of Maxwell's theory in four spacetime dimensions. 
The non-singular nature of the solution motivated Melvin to dub his solution a magnetic 'geon'. 

For the generalised Melvin solution, the existence of singularities will obviously depend on the 
choice of Lagrangian. From the Einstein equations (|2.5D one reads for the Ricci scalar 



R = 8 



L(s) - s^l 
as 



(4.8) 

H2 



2A(p)4 



so that if the energy density or its gradient blows up so does the curvature. For the Euler-Heisenberg 
family this curvature invariant is bounded everywhere, but for the Born-Infeld case it diverges at 
radial distance 

A(p) 4 = P 2 H 2 . (4.9) 

Therefore, if the constant H in the solution is such that H 2 > /3 -2 , then the solution is singular. 
The value H 2 = f3~ 2 is the maximal magnetic field for a purely magnetic configuration in Born- 
Infeld theory. Notice that Born-Infeld theory has a maximal value for the electric field E and 
for the magnetic field H (but not for the inductions). Thus, the free parameter describing the 
solution is constrained if we require absence of curvature singularities. In general, we expect that if 
the non-linear theory of electrodynamics under consideration has a critical field the curvature will 
diverge at that critical field. 

Independently of the existence of curvature singularities, the metric ( |3.1| ) will be singular if 
there exists a radial distance po such that 

pi + ^h(po) = 0. (4.10) 

This is certainly the case for Euler-Heisenberg and Born-Infeld, since h(p) is negative everywhere 
for these two theories. Consider the p-cj) sections of the metric ( |3.1| ) 

o 2 

ds 2 = ^—dp 2 + g H d<t> 2 . (4.11) 
9<t><t> 

In the neighbourhood of p = po, 

ds 2 -—^ d p 2 + g H {p-p )d4> 2 , (4.12) 

9<M>{P ~ Po) 

where = (dg^/dp) at p = po, which we assume to be non-zero. Now introduce the new 
coordinate R by 

R 2 

(P ~ Po) = 4 500 — , (4.13) 
Po 
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to get the metric 

ds 2 ~dR 2 + R 2 d[^^<py (4.14) 
Hence, if we identify <f> with period 

- np0 (4.15) 



the singularity at p = po is just the usual singularity at the origin of flat space in polar coordinates. 
Any other period for <f> leads to a conical singularity. In any case, the point p = po should be faced 



as the location of the axis of symmetry, where po is the first solution of (4.10) when coming from 
infinity. The region < p < p will be discussed below. 

Asymptotically, our generalized Melvin Universe will approach the usual Melvin solution, which 
is clearly a consequence of the requirement that the theory of non-linear electrodynamics under 
consideration should approach Maxwell's theory for weak enough fields. The components of the 
Riemann tensor in an orthonormal frame vanish as 1/p 6 . In this sense, asymptotic Melvin space- 
times are asymptotically flat (away from the axis of symmetry). The components of the Weyl tensor 
in an orthonormal frame, C a bcd , also vanish with the same power of p. The fact that this is true for 
all non-vanishing components indicates that there might exists a conformal compactification of the 
spacetime, which however is not know. If it exists the conformal spacetime should have a light-like 
conformal boundary. 



We also note that the isometry group of the solution (34) is the same as the one of ( |3.7D . There 
is a U(l) x ISO(l, 1) isometry corresponding to a circle action along the <j) direction, boosts along 
the t-z direction and translations along the t and z directions. 

4.2 Test motions and supersymmetry 

In the Newtonian limit, the generalised Melvin Universe has a gravitational potential equal to the 
usual Melvin solution, which is of the form 

tt2 tt4 

V = fLf + fL p \ (4.16) 
4 h 32 H \ ) 

This potential well will tend to confine the particles in a cylinder of finite radius. In this respect 
the behaviour is similar to that of Anti-de-Sitter spacetime, despite the fact that AdS has constant 
curvature in contrast to the asymptotically vanishing curvature of Melvin discussed in the previous 
subsection. For the full relativistic analysis consider the geodesic equations, where A is an affine 
parameter along the trajectory: 



dp\ p 2 + 4H- 2 h(p) 



dX J p 2 A(p) 



2 



1 /77!2 2\ 2 

(h — p ) — m 



Hp) 



(4.17) 



lA(p) 2 dt E dz p 



dX p 2 + AH- 2 h{ P y dX A(/o) 2 ' dX A(p) 2 ' 

The conserved quantities E,p,l are associated with the Killing vector fields d/dt,d/dz and d/d(j) 
respectively. The 'mass' m is as usual associated with the metric Killing tensor. The t, z coordinates 
are the proper time and proper radial distance for an observer at the origin. As seen by this observer 
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clocks tick faster as the radial distance increases. This is exactly the opposite of the redshift 
witnessed by the observer at infinity for events close to the black hole horizon in a Schwarzchild 
spacetime. 

From the radial equation we can see that at large distances the repulsive terms dominate over 
the attractive (energy) term. Hence, no particles may enter the large enough p region, whatever 
values the particle's quantum numbers take. The solely exception is a photon with no angular 
momentum not moving entirely in the z direction, which is not confined to a cylinder of some 
radius. In this way we see that the Newtonian intuition is mostly confirmed. This resembles the 



Melvin case [15] due to the similar asymptotics. 

The behaviour for small p is quite distinct from the usual Melvin solution. From the radial 
equation we can see that no particle can cross the p = po cylinder. But geodesies may exist on 
both sides of this cylinder. As an example consider the Euler-Heisenberg case. Equation ( |4.10| ) 
only has one solution. Its left hand side is negative in < p < po and positive for p > p$. Bound 
states with E 2 < m 2 can therefore exist within the p = po cylinder. In a sense this is analogous to 
the r < region in the Kerr metric. Although p = pq should be seen as the location of the axis 
of symmetry (like the r = in Kerr), there is a region beyond it for which the metric is different 
and therefore taking it seriously is not double counting (like the negative p region for flat space in 
polar coordinates). Hence the spacetime can have two geodesically complete disjoint regions. 

To close this section we remark that on a curious effect for the Einstein-Maxwell Melvin Uni- 
verse: at any radial distance one may find a 'no-force' configuration, i.e. a static test particle. The 
charged geodesies are 



d\ J A(p) 4 ( E P ) A(p) 2 p 2 V BA(p) J 



(4.18) 

d0 A{pf f _ 2q \ 
dX p 2 V BA(p))> 

where q is the particle's charge. The t and z equations are the same as in ( pLTTj ). Notice that the 
charged geodesies for the generalised Melvin Universe are more subtle due to the coupling to the 



Maxwell potential. Choosing I, q as to cancel the right hand side of the 4> equation of motion (4.1 



for a given p, putting p = to cancel z motion Q4.17 ) and choosing E 2 = A(p) 2 m 2 , the particle 



becomes static. Obviously this cannot happen for photons {m = 0). 

The possibility of a no-force configuration, together with the fact that the Reissner-Nordstrom 
family includes a supersymmetric solution raises a question about the supersymmetry of the Melvin 
solution within = 2, D = 4 supergravity. The Killing spinor equation is, in the notation of forms, 

d + \u ah T ah - \F ab T ah T^ e = 0. (4.19) 

We use for the flat metric 7] a b = diag(—l, 1, p 2 , 1), i.e., flat space in cylindrical coordinates. Then, 
one easily checks that the Killing spinor equation implies the conditions: 

e = A5e , (l-^vAe = 0, (4.20) 
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where eo is a constant spinor. Keeping in mind that {r^,r^} = 2p~ 2 one concludes that the 
right hand side equation has no solution all over the spacetime. However for any value of B 
( [4.20D is obeyed on the hypersurface B 2 p 2 /4 = 1. This hypersurface 'preserves half of the vacuum 
supersymmetries'. Note that if we require the classical expression for circular orbits of a charged 
particle in a constant magnetic field, p 2 = l/(qB), where B is the magnetic field measured at radial 
coordinate p, the no-force condition I = 2qj (BA(p)) singles out this 'supersymmetric' hypersurface. 
In a subtle way, therefore, the no-force condition is still associated with unbroken supersymmetry 
of the solution. The phenomenon uncovered here resembles what are called hypersurface twistors 



4.3 Dirac Quantization of the Magnetic Flux 



Inverting the constitutive relations (2^3), one finds 



F llv =L s P llv + L t *P liV . (4.21) 

We will assume that Lt = for a purely magnetic configuration, which is true for most theories of 
interest (including Born-Infeld and Euler-Heisenberg). Then 

F„ = Ls*. - ^ = -^j( 1 -"W + ^)* ^ 2 > 

where F = dA. The subscript 00 indicates this is a good gauge at infinity. This is not the case at 
po as defined by ( 4. 10| ) . A good gauge at po is Aq = Aoo + d\ with 



2 / 1 dh, A / . 

For fermions minimally coupled to the Maxwell field, this gauge transformation induces a change 
of phase in the wave function required by gauge invariance, of the form fy(x) — > exp (iex)^(x), 
where e is the electron charge. To be singled valued we have the following quantization condition 
for the constant H 

1 dh . A n ,. 
1 + — -^-r (po) =tt- (4-24) 



2ttH V 2p H 2 dp K ' ' J 2e 

The left hand side of this equation turns out to be the magnetic flux (of the magnetic induction 
B) through a z = const, t = const surface, 

®b = ^T F^dS^ = — / Frfdfxty. (4.25) 
71 " ./ ^ J 0=0 Jp= P0 



Thus, the flux for the magnetic induction is quantized, and one reads off the usual Dirac quantization 
condition, 

where we have reinserted fundamental constants. In the Einstein-Maxwell case, the analogous 
condition to ( 4.24] ) takes the form 2e = Bn, and one gets the universal quantization condition for 



the flux from the fact that &b = B .In the generalised case, the metric parameter H obeys the 
more involved constraint ( 4.24f) but the universal flux quantization still holds. 
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5 Travelling Waves in the Magnetic Universe 



A spacetime solution to Einstein-Maxwell theory describing waves of arbitrary profile travelling 
in a magnetic universe was found in [||. The solution is described by the fields {g^v^F^v) and is 
obtained by using a Kerr-Schild ansatz: 

9[iv = 9iiu + A(p)" 2 ^(n, p, (fyk^ky. (5.1) 



The fields ^g ilu -,F ilu ) are the ones in the usual Melvin solution ( |3.7| ), and the vector k = k^d^ = 
d/dv is null and Killing for the Melvin geometry. We have introduced light-cone coordinates 
(u, v) = {z rfc t) /\/2- The Einstein-Maxwell equations are then still solved if the function ^(u, p, <p) 
is harmonic in the usual Melvin spacetime, i.e. D^D^^(u, p,(p) = 0, where is the covariant 
derivative for <7„„. The harmonic condition does not constrain the u dependence (hence the arbitrary 
profile described by an arbitrary function f(u)), whereas the </> and p dependence might be solved 
by separating variables. The geometry obtained reads 

o 2 d6 2 

g^dx^dx" = A(p) 2 (2dudv + dp 2 + f(u) cos m{<j> - (f> )P (j>) du 2 ) + fj—^ . (5.2) 

The function P{p) obeys a second order ordinary differential equation which can be computed 
to arbitrary accuracy by an iteration method described in ||, whereas 4>o and m are constants. 
The waves are travelling along the k direction which lies on the Einstein light cone for the Melvin 
Universe. 

It is straightforward to generalise this solution to the non-linear case. Using again the ansatz 
( phi] ) for the metric, where now g^ v is given by ( |Q| ) one can check that (g^, -P^jy) is a solution 
to (|2.5|), where P^ u is still given by (|3~l|), as long as the function ^>{u,p,<j)) is harmonic in the 
geometry of the generalized Melvin Universe. By separating variables we find that *$>(u,p,(f)) = 
f(u)R(p) cosm((J) — 4>q) where R(p) obeys: 



d_ 

dp 



y+jg^M MMN _ ™w m _ . (5 . 3) 

v p dp J p 2 + AH~ 2 h{p) K ' K ' 

For large p this approaches the solution discussed in |J. The geometry describing waves in the 
generalized Melvin Universe is then 

ds 2 = A(p) 2 (2dudv + f(u) cos m(0 - ^o)^ (p) du 2 ) + p2 fw-*h(p) + ~ + T(p) 2 ~ {5A) 
To conclude this section, we remark that the metric describing the propagation of fluctuations 



in non-linear electrodynamics, the Boillat metric, is not conformal to the Einstein metric [12]. The 
fact that the waves in the generalized Melvin Universe travel along a null direction of the metric 
( |3.1| ) is because this is the principal null direction of F^, and therefore along this direction the 
Einstein and Boillat light cones coincide. This would not be the case anymore if one could find 
solutions to (p.5|) of the Kerr-Schild type, where k would have non-zero (j) or p components. 



6 Electric-Magnetic Duality 
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6.1 Einstein-Maxwell 



The limiting procedure we described in section 3 also works if we start from the dyonic Reissner- 
Nordstrom black hole. In the usual stereographic coordinates this solution is 

4r 2 



ds 2 



-Vdi 1 + V~ l dr 2 + 



(1 + CC) 



F 



Q 



dtAdr-P sin 6d6 A , 



At At T/ 1 2M j_ Q 2 + P2 

; d(d(, V = 1 + = . 



Q , , „ d(Ad( 

-^dt Adr- 2iP—^ 

r 2 1 + CC 2 



(6.1) 



Performing the method of section 3.1, taking also the limit P — > oo with PA 2 = fixed = p, and 
performing a final set of coordinate transformations 



t 



E 2 + B 2 ' 



T 

x = l r 



z 

r r = 1 + 



— - A — p = Hp), 



3.2) 



q = iB, p = iE, m = —(E + B )/2, 
we then end up with the dyonic version of the Melvin Universe 

ds 2 = p 2 A(pr 2 d<p 2 + A(p) 2 [-dT 2 + dZ 2 + dp 2 } 

F = EdT AdZ + BpA(p)- 2 dp A d<p. 



(6.3) 



Of course (|6.1|) and (6.3) can be obtained by SO(2) duality rotations of Einstein-Maxwell theory in 



four dimensions from (3.3) and (3.7) respectively. Therefore the limiting procedure comutes with 
electric-magnetic duality. 



6.2 Non-linear electrodynamics 

For theories obeying the requirement F^-kF^ 11 = P^-kP^, there is also a electro-magnetic duality 
rotation that leaves invariant the equations of motion and energy- momentum tensor [11]. One 
therefore obtains new solutions to (p.5|). Born-Infeld is one such theory. We can thus obtain the 
dyonic Melvin Universe for Born-Infeld: 



ds 2 = p2 + ^ Hp) d6 2 + 



p 2 Hp) 2 



MpY 



P 2 + tjt+siHp) 



dp 2 + A(p) 2 (-dT 2 + dZ 2 ), 



(6.4) 



P = - 



Hp 

Hp) 2 



with 

Kp) 



^-J (1 +n)!(l + 4n) 

n=l v ' y ' 



[H 2 + D 2 )(3 : 
2A(p) 4 



A dp + DdT A dZ, 



[](2m-l) 

,m=l 



. , x D 2 + H 2 r> 

A(p) = l + ; P 2 . (6.5) 



For theories without duality invariance one can obtain the Melvin electric universe starting with 
the magnetically charged Reissner- Nordstrom solution in the theory (see [|llj) and applying the 
limiting procedure. 
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7 Discussion 



The limiting procedure we have described is very simple both mathematically and conceptually. 
However, notice it changes completely the causal structure and singularity structure of the solution: 
in the linear case, i.e. Maxwell's theory, we start with a singular solution with a null surface and 
we end up with a completely non-singular solution with an everywhere timelike Killing vector field. 
With the wisdom of hindsight the procedure looks quite natural; but it is far from obvious that a 
charged black hole and a magnetic flux tube are so simply related. 

In our generalisation of the procedure to the case of non-linear electrodynamics, taking the limit 
commutes with solving the equations of motion. However, one should not take this for granted. 
For instance, applying this procedure to a Reissner-Nordstrom-Anti de Sitter (or de Sitter) black 
hole spacetime , one ends up with 

F = BpA(p)~ 2 dp A d(f>, 

where is the cosmological constant. In this case the equations of motion are only solved if 
B 2 = 1. The procedure might, however, be useful in finding other solutions of the flux tube type. 
It is perhaps worth exploring applications to black p-branes of string theory. 
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